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ABSTRACT 

Recent observations show that hypernovae may deposit some fraction of their kinetic energy in mildly 
relativistic ejecta. In the dissipation process of such ejecta in a stellar wind, cosmic ray protons can be 
accelerated up to ~ 10 19 eV. We discuss the TeV to MeV gamma-ray and the X-ray photon signatures 
of cosmic rays accelerated in hypernovae. Secondary X-ray photons, emitted by electron-positron pairs 
produced via cascade processes due to high-energy protons, are the most promising targets for X-ray 
telescopes. Synchrotron photons emitted by protons can appear in the GeV band, requiring nearby 
(< 40 Mpc) hypernovae for detection with GLAST. In addition, air Cherenkov telescopes may be able 
to detect regenerated TeV photons emitted by electron-positron pairs generated by CMB attenuation 
of 7T° decay photons. 

Subject headings: supernovae: general — gamma rays: theory — X-rays: general — radiation mech- 
anisms: nonthermal — cosmic rays 



1. INTRODUCTION 

Hypernovae (HNe) are a peculiar type of supernova 
with ejecta velocities and apparent isotropic-equivalent 
ejecta energies which are larger than usual, and with 
indications of anisotropy (Nomoto, et al, 2008). Some 
of them are associated with long gamma-ray bursts 
(e.g. Woosley and Bloom, 2006; a recent example be- 
ing GRB060218/SN2006aj, e.g. Campana et al, 2006), 
while others appear not to be. Cosmic rays up to en- 
ergies -E<,10 17 eV are thought to be accelerated in rela- 
tively normal supernova remnants (Hillas, 2005). More 
recently, it has been suggested that HNe may accelerate 
cos mic ray protons or nuclei up to energies E <Ax 10 18 Z 
eV (IWang et al. ll2007HBudnik et aT1l2007l:rWang et al. I 
l2008f h The smoking gun proof for cosmic ray accelera- 
tion in supernova remnants would be the observation of 
very high energy (>, TeV) neutrinos (e.g. Kistler and 
Beacom, 2006), which will require completion of cubic 
kilometer detectors, or else the observation of secondary 
photons arising from pions, which remains inconclusive 
(e.g. Gabici and Aharonian, 2007; Katz and Waxman, 
2008; etc.). Nonetheless, continued gamma-ray observa- 
tion with air Cherenkov telescopes (ACTs) in the TeV 
range, and with GLAST and AGILE in the GeV range 
may provide the best immediate hopes for resolving this 
question. In the present paper, we consider the same 
question in relation to HNe, and address the question 
of the photons signatures from secondaries arising from 
cosmic ray acceleration in typical HNe. 

In |j2l we describe our model NHs and consider the 
baryonic and photonic environment in which the explo- 
sion occurs, as well as its effect on the cosmic rays accel- 
erated in the ejecta. In fj3] we discuss the Monte Carlo 
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simulations performed on these models, and present the 
results for the photon signatures arising from various sec- 
ondary components. In [|4] we discuss the detectability 
prospects for these signals, compared to the sensitivity 
of ACTs and GLAST, and summarize our results and 
conclusions. 

2. MODEL DESCRIPTION 

HNe, especially the Ic types associated with GRB but 
also some of the unassociated ones, are thought to be due 
to WR progenitors (Nomoto et al, 2006), and as such are 
expected to have had a strong stellar wind phase prior to 
the explosion. The model of a HN ejecta expanding in a 
stellar wind used in this paper is based on the model of 
IWang et al. I (p007n . Thus, we consider a stellar wind en- 
vironment around the progenitor which is characteristic 
of WR stars. Assuming a mass loss rate M and a wind 
velocity v w , the density profile of the wind is written as 
p(r) = 5 x 10 11 A»r _2 g cm -1 , where A* (oc M/v w ) = 1 
corresponds to M = lO -5 M0yr -1 and v w = 10 3 km s _1 
(|Wang et al. Il2008h . 

The outer envelope regions (the ejecta ) of th e explod- 
ing HN, as shown by ISoderberg et al. I (|2006n . have a 
kinetic energy distribution oc (T(3)~ a . In this paper we 
assume E k = 10 52 (r/3/0.1)" 2 erg, where the velocity of 
the bulk of the ejecta ranges from Tf3 ~ 0.1 up to semi- 
relativistic values {Tj3 ~ 1), where (3 = v/c and T are 
the ejecta normalized velocity and bulk Lorentz factor, 
respectively. 

For non-relativistic ejecta (r/3 < 0.5, £ k > 4 x 10 50 
erg), the free expansion phase before deceleration sets 
in lasts for 440(r/3/0.5)~ 5 ^4~ 1 days. Therefore, non- 
relativistic ejecta cannot dissipate their kinetic energy 
within the 10-20 days typical timescale of the UV-optical 
photon radiation from HNe. Thus after the optical emis- 
sion from HNe has declined, we may not expect sec- 
ondary photons originating from ^-interactions, even if 
a sufficient amount of high-energy protons are produced. 

On the other hand, the bulk of the kinetic energy of the 
mildly relativistic ejecta are dissipated at a radius i?d ~ 
lO^r/Vl.O)- 1 ^- 1 cm. Since /3 ~ 1 for such ejecta, all 
the mildly relativistic ejecta dissipate their kinetic energy 
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within 10 days, as long as A* > 1. Therefore, hereafter 
we consider only the most energetic component T/3 = 1 
(Ek — 10 50 erg) of the mildly relativistic ejecta. 

For the ejecta with Tf3 — 1, the magnetic field at the 
dissipation radius i?d may be estimated as B 2 /8tt = 
€Bp{Rd)c 2 , where £b = O.le^-i is the fraction of the 
equipartition value of the magnetic field. Thus, we ob- 
tain B ~ $.4e)£_ x A i J 2 G. 

The maximum energy of accelerated protons for the 
mildly relativistic ejecta may be written as £ max — eBR. 
With the above values of B and i?d, we obtain e max ~ 
10 19 e 1 J 2 _ 1 Al^ 2 eV, for the most energetic ejecta dissi- 
pated within 10 days. Adopting this value, the en- 
ergy distribution of accelerated protons is assumed to be 
N p (e p ) oc e~ 2 exp (— e p /e max ). The amount of protons 
is normalized by an efficiency factor 1/6 for the conver- 
sion of dissipated ki netic energy ( 10 50 erg) into energy of 
accelerated protons (|Hillas 1 120051 ). 

The usually observed low energy (~ 1 eV) photons 
from HNe are attributed largely to radioactive decay in 
the non-relativistic (Tf3 ~ 0.1) ejecta. After an initial rise 
and just after the dissipation of the ejecta (~ a few days), 
the changes in the optical luminosity of the HN are not 
very d rastic (1-2 magnitudes), as shown in iPian et al. I 
(2006). Therefore, a constant luminosity for a few tens 
of days may be a reas onable approximation. As shown in 
iMazzali et al~l (j2006l ). the peak flux is at 4000-5500 A, 
which is considered to be due to broad emission lines of 
several types of metals. Here we mimic this photon field 
by a thermal photon field with a temperature T = 1 eV 
and luminosity Lhn = 10 erg s _1 , which is compatible 
with the bolome tric luminosity of SN 1998bw around its 
peak brightness (|Galama et al. IIT998D . 

Another source of low energy photons is synchrotron 
radiation from electrons accelerated in the ejecta. For 
mildly relativistic shocks, the typical Lorentz factor of 
accelerated electrons is less than m p /m c ~ 10 3 , which 
implies that the typical energy of synchrotron photons 
is e 7 < 10 -2 eV for B = 3 G. The cooling timescale of 
such electrons is ~ 1 day, which is much shorter than the 
timcscales we consider. Injection of accelerated electrons 
may continue for a few tens of days as indicated by radio 
observations. The most luminou s radio afterglow i n HNe 
observed ever is in SN 2003dh (|Berger et al. ||2003D of • 
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10 41 erg s _1 , which may slightly enhance pion production 
efficiency for the highest energy protons (> 10 18 eV). 
Here, for simplicity, we neglect the photon emission from 
accelerated electrons. 

In order to calculate the pion production and the sub- 
sequent cascade processes, we carry out Monte Carlo sim- 
ulations using a code developed for and dis cussed in a se- 
ries of studies of gamma-ray burst physics (Asano Il2005t 
lAsano fc Nagataki II20M IA sano & Inoue 1120071 ). We use 
a simple one-zone approximation, and follow during a fi- 
nite time the physical processes of pion production, pion 
decay, muon decay, and 77 electron-positron pair cre- 
ation, as well as the usual radiation processes of syn- 
chrotron (SY) and inverse Compton (IC) emission from 
protons, pions, muons, electrons, and positrons. We as- 
sume that the accelerated protons are injected promptly 
just after the energy dissipation occurs at r = R^. 
Since the bulk motion of the post-shock region is non- 
relativistic, we neglect the expansion of the ejecta after 
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First, we consider a HN exploding in a standard wind 
with A* = 1. As shown in Fig. [TJ even after the acceler- 
ated protons have been irradiated by the HN photons for 
Rd/c ~ 4 days, most of the protons still retain a large 
fraction of their energy. The pion-production timescale 
i P7 is inversely proportional to the HN photon density 
cx LT- X R~ 2 . Following IWang et al. I (|2007l) . the frac- 
tion of the energy lost by protons to pions is estimated 
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Fig. 1. — Spectra of the initial protons, and of the surviving 
protons/neutrons after the timescale Rd/c = 3.3 X 10 5 s ~ 4 days, 
for a hypernova explosion in a wind characterized by the mass loss 
rate A* = 1. The spectra of secondary particles and photons cre- 
ated during this time scale are also plotted, including pions, proton 
synchrotron photons (labeled p-SY), and proton inverse Compton 
photons (labeled p-IC), not all of which can be directely observed. 
The fluence is normalized assuming a distance D = 100 Mpc. The 
hypernova soft photon spectrum is plotted with thin dashed lines. 



As shown in Fig. [2l if one neglects the 77-absorption 
effects, a spectral bump of SY emission appears around 
the TeV-10 TeV energy range. This bump is due to SY 

photons of 15e 1 J 2 _ 1 A^ 2 TeV from positrons and electrons 
from muon decay. Since the SY photon energy scales as 
the square of the energy of the emitters, the photon bump 
energy range spreads over 2-3 orders of magnitude. The 
inclusion of 77 effects leads to these photons being ab- 
sorbed through interaction with the HN soft photon field, 
creating electron-positron pairs of 1-10 e 1 J 2 _ 1 A^ 2 TeV. 
The secondary pairs emit SY photons over a wide en- 
ergy range (which should be twice as wide as that of the 
TeV bump), from a few keV to a few hundreds of keV. 
These pairs cool so promptly that the result is a sim- 
ple power-law spectrum due to cooled pairs below 100 
eV. The secondary photon emission is expected to last 
as long as the HN emits optical photons, even though 
the cooling time of the electrons that emit X-ray pho- 
tons is only about 1 minute. These photons will be ob- 
served with various present-day X-ray telescopes for 100 
ks integration, as shown in Fig. [2l unless the X-ray af- 
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Fig. 2. — Spectral fluencc of combined secondary and tertiary 
photons from the hypernova of Fig. [T] (thick solid line), at the dis- 
tance D = 100 Mpc for A* = 1, integrated over R d /c = 3.3 X 10 s 
s. Thin curves denote electron/positron synchrotron (labeled e- 
SY), proton synchrotron (labeled p-SY), and 7r°-decay components 
emitted inside the source. The resulting photon spectrum escaping 
from the source, taking into account internal 77 absorption by the 
hypernova soft photon field, is shown by heavier lines. The sensi- 
tivity curves for representative instruments are plotted with dotted 
lines. The curves for XMM-Newton and Suzaku are normalized by 
100 ks integration time. 

Around the GeV region, proton SY emission yields 
a bump in the photon spectrum, as seen in Fig. [2l 
The typical energy of proton SY photons is 2Ae 3 J 2 _ 1 A^ 2 
GeV, and the ratio t pl /t PiSyn (where t P}Syn = SY cooling 

timescale) is ~ 0.2(A»es,-i) 3 ^ 2 , which roughly agrees 
with the obtained energy fraction of proton SY to pi- 
ons. For this distance the fluence is much lower than 
the GLAST detection limit. To detect 100 MeV photons 
with GLAST would require a HN at < 6 Mpc. 

Another notable feature of Fig. [2 around energies 
~ 10 16 eV, is the prominent presence of photons from 
7r°-decay. Those photons escape without being absorbed 
by the ~ 1 eV thermal photon field assumed for the HN. 
However, the mean free path of these 10 16 eV photons 
against 77 absorption by co smic microwave backgro und 
(CMB) photons is ~ 10 kpc (jAharonian et al. |[200l ). so 
that we cannot expect to detect such photons directly. 
The secondary electron-positron pairs generated by at- 
tenuation are very energetic, and are inverse-Compton 
scattered by CMB pho t ons, e.g. as discussed for GRBs 
(jRazzaque et al. I l2004t iMurase et al. I I2007L and refer- 
ences therein). These boosted photons can pair-produce 
again, and the process repeats itself until the energy of 
the degraded photons is in the 1-10 TeV range. The 
mean free path of these regenerated 1-10 TeV photons is 
longer than 100 Mpc, and they can reach the Earth. As 
long as the intergalactic magnetic field is weak enough, 
the delay time of TeV photons emitted by ~ 100 TeV 
electrons/positrons is negligible in compariso n with the 
timescale a few days ()Razzaque et al. I [2004V We omit 
plotting the spectrum of the regenerated photons in Fig. 
O since from D = 100 Mpc it will be hard to detect them; 
however, if a HN occurs in the Virgo Cluster (D ~ 20 
Mpc) , there would be a chance to detect these secondary 



TeV photons (see below). 

Next, we consider a HN occurring in a denser wind 
with A* — 5. This value is compatible with cur- 
rently available data on wind mass losses suggesting 
M < few x 10~ 4 Mfp/yr, which refer to st ars well be- 
fore any explosion (|Mevnet &: Maeder 120071 ). Physically, 
even larger values may be plausible, since one expects 
the the mass loss to increase considerably as the evo- 
lution of the core rapidly approaches the final collapse, 
with a rapid increase in the luminosity and the envelope 
expansion rate. 

In this case of A* = 5, the ejecta will stall at Rd ~ 2 x 
10 15 cm within ~ 1 day, so that the non-relativistic ejecta 
of Tf3 = 0.1 can catch up with the decelerating ejecta 
about 10 days later. However, at least until the non- 
relativistic ejecta has caught up, the secondary photons 
are largely observable. 

Basically, the cooling time-scale is shorter than the in- 
tegration time scale (~ 4 days) assumed here, which re- 
sults in a bumpy non-power law energy distribution of 
the final protons (Fig. [3]). The highest energy protons 
cool via SY (t P7 /t p ^ syn ~ 1), while protons of < 10 18 
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Fig. 3. — The primary and secondary particle and photon spec- 
tra, similar to Fig. [T] for a hypernova at the same distance D = 100 
Mpc but in a denser wind environemnt of A* = 5. 

In this case the proton SY emission becomes promi- 
nent, since the cooling timescale cx £ B 3 L\ A* 7 ^ 2 is shorter. 
The secondary photon flux (see fig. 0} is, as expected, 
larger than in the lower density wind (fig. ^ case. 
However, even this higher density wind case gives, from 
D=100 Mpc, an insufficient flux to be detectable with 
GLAST (see Fig. [4]). However, photons from similar 
HNe within 40 Mpc would be detectable by GLAST. In 
addition, the regenerated TeV photons are promising tar- 
gets of ACTs. In Fig. [4] we plot the regenerated photon 
sp ectrum obtained by t he same numerical simulation as 
in IMurase et al. I ((2007) . They are well above the detec- 
tion limit of present-day ACTs even for D = 100 Mpc. 

Secondary X-ray photons are emitted by electron- 
positron pairs originating from photons due to both pro- 
ton SY and SY photons from muon-decay positrons. This 
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Fig. 4. — The secondary and tertiary photon spectra, similar to 
Fig. [2] but for the hypernova of Fig. [3] at D = 100 Mpc in a wind 
of At, = 5. In addition, the muon synchrotron component (labeled 
/i-SY), whose typical energy is determined by a balance between 
the cooling timescale and lifetime (sec e.g. Asano 2005), is also 
plotted. The dash-dotted line is the regenerated photon spectrum 
due to CMB attenuation of ~ 10 16 eV photons. 

is because the energy of the proton SY photons shifts 
higher oc e 3 J 2 _ 1 A^ 2 . Since the energy range of the ab- 
sorbed proton SY photons (~ 10 11 eV) is narrow, the 
secondary pairs produce a characteristic peak in the flu- 

1 /2 3 /2 

ence spectrum around 0.4e B ' _ l AJ keV. These photons 
are indirect evidence of proton SY, and can be easily 
detected with present-day instruments. 

4. SUMMARY AND DISCUSSION 

We have shown that secondary gamma and X-rays, 
correlated with the initial thermal optical emission of 
hypernovae in the first ~ 10 days, can provide evidence 
for proton acceleration, as well as provide a diagnostic 
for amplification of magnetic field in the blast wave, and 
for the mass loss rate in the progenitor stellar wind prior 
to the explosion. There are three main spectral compo- 
nents of secondary photons: 1) X-ray photons emitted 
by electron-positron pairs originating from 77 interac- 
tions initiated by synchrotron photons from muon-decay 
positrons or protons, 2) synchrotron photons emitted 
by protons in the GeV band, and 3) regenerated TeV 
photons emitted by electron-positron pairs generated by 
CMB attenuation of 7r° decay photons around 10 16 eV. 
The X-ray photons are the most promising targets, so 
that follow-up observations of HNe with X-ray telescopes 
are indispensable to find evidence of proton acceleration. 
Soft SY photons from accelerated electrons may also en- 
hance the electromagnetic cascades by interacting with 
photons from pion-decay. The interesting 7r°-decay pho- 
ton signature (component 3) is also an interesting can- 
didate for detection with ACTs in the dense wind case 
(A* = 5). 

If a HN occurs in our Galaxy at a distance of 10 kpc 
(the rate for which should be ,<10 _3 — 10~ 4 yr _1 ), our 



results indicate an expected flux 10~ 7 -10 -5 erg cm -2 
s -1 at 10 GeV for A* = 1-5, due to proton SY and/or 
secondary leptons, detectable by GLAST. By compari- 
son, the most luminous "norm al" SNRs observed with 
EGRET (jEsposito et al. Ifl996h have fluxes of - 10~ 10 
erg cm -2 s _1 at 10 GeV. TeV photon detections are not 
expected to be detectable from a Galactic HN, since the 
photon regeneration process mean free path is too long 
to be effective here. Our simulations show also that the 
secondary ^ TeV neutrinos from the cascades in a Galac- 
tic HN have a spectral peak at 10 16 eV with a flux of 
10~ 4 -10 -3 erg cm -2 , well above the detection limit of 
IceCube. Thus, one would expect to detect continuous 
TeV neutrino emission for a few days from such Galactic 
HNe. 

For HNe at distances D ~ 100 Mpc, the X-rays (com- 
ponent 1) will be easily detectable by XMM, the sub- TeV 
radiation is marginally detectable near the low energy 
threshold by MAGIC and similar ACTs, and the GeV 
photons from a proton SY (component 2) arc difficult to 
detect with GLAST. However, for A* = 5, HNe in the 
Virgo Cluster (D ~ 20 Mpc) would be easily detectable 
at GeV energies by GLAST, as would also the regener- 
ated TeV photons. Such detections would provide con- 
straints on HN models. E.g. the duration of the proton 
SY emission or the X-ray spectral peak at 1 keV gives 
an estimate of the survival timescale of magnetic fields 
amplified by the non-linear MHD turbulance excited by 
cosmic rays (Bel l fc Lucek1l200l . Since the maximum 
proton energy is not so sensitive to eb, the spectral com- 
ponent 1 due to the cascades from pion production will 
not change drastically, even for eb <C 0.1, although the 
direct proton SY (component 2) can become negligible. 
In such cases, the disappearance of the X-ray spectral 
peak shown in the A* = 5 case would be a diagnostic for 
such low eb values. 

The intensity of the cascades depends on the wind den- 
sity, providing a diagnostic for the progenitor mass loss 
rate. For example, with A* = 0.2 and other parameters 
as in Fig. [2] the fluxes are undetectable even by X-ray 
instruments, unless the source is extremely near. On the 
other hand, larger values of A* ~ 10 (e.g. as suggested by 
Campana et al (2006) for SN2006aj) would give higher 
fluxes than those of Fig. 0] enhancing the probability 
of detection at D ~ 100 Mpc at TeV, GeV, and X-ray 
energies. 
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